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£^ ' Abstract. - We present a new dispersion relation for photons that are nonlinearly interacting 

with a radiation gas of arbitrary intensity due to photon-photon scattering. It is found that the 
■ photon phase velocity decreases with increasing radiation intensity, it and attains a minimum 

value in the limit of super-intense fields. By using Hamilton's ray equations, a self-consistent 
kinetic theory for interacting photons is formulated. The interaction between an electromag- 
netic pulse and the radiation gas is shown to produce pulse self-compression and nonlinear 
saturation. Implications of our new results are discussed. 

a,: 

qj i In classical electrodynamics photons are indifferent to each other in vacuum, while in 

• quantum electrodynamics (QED) photons can interact via virtual electron-positron pairs, 
giving rise to vacuum photon-photon scattering [1]. This is commonly expressed using a 
series expansion (in the field strength) of the Heisenberg-Euler Lagrangian, yielding nonlinear 
corrections to Maxwell's vacuum equations. These corrections give rise to single particle effects 
(some of which are highly speculative) , such as closed photon paths [2] , vacuum birefringence 
[3-5], photon splitting [4] and lensing effects in a strong magnetic field [6,7], and coherent 
effects, such as the self- focusing of beams [8] or the formation of light bullets [9]. Moreover, 
as an example of a collective effect, the QED photon self-energy will change the refractive 
index of a radiation gas [10-12], which can result in, e.g., Cherenkov emissions in a radiation 
gas [13]. 

Recently, it has been shown that nonlinear vacuum corrections cause photonic collapse in 
a radiation gas [14]. The phenomena of photonic collapse in two space-dimension has been 
confirmed by computer simulation studies [15]. As indicated in Ref. [15], the increase in 
intensity will show no upper bound as the collapse has begun, and will eventually surpass the 
Schwinger field ~ 10 16 V/cm 2 . Even as dispersive corrections are added [16-18], with resulting 
pulse splitting, the field intensities within the collapse region may still reach values well above 
the ones allowed by the approximate Heisenberg-Euler formalism. Thus, this raises the critical 
question of how far the predictions of the approximate Heisenberg-Euler Lagrangian can be 
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trusted, i.e. whether the collapse of electromagnetic fields is a "true" effect or an artifact of 
the approximation used. 

The investigation of the occurrence of singularities is important not only because of its 
theoretical interest, but also by the fact that Nature may offer the environment for such 
events to take place. Within astrophysical and cosmological settings, such as neutron stars 
and magnetars [19], the conditions are met for the ensue of collapse of disturbances in the 
photon distribution. These conditions may also soon be met in laboratory environments. 
Today's lasers can produce 10 21 — 10 22 W/cm 2 [20], and the next generation laser-plasma 
systems are expected to reach 10 25 W/cm 2 . [21-23]. This has been used to suggest that 
field strengths close to the Schwinger critical value could be reached [24]. Thus, here the 
nonlinear effects introduced by photon-photon scattering will be significant, and the weak 
field approximation no longer holds. 

In this Letter, we investigate a photon gas of arbitrary intensity at the one-loop level. 
We begin with a discussion of the dispersion relation for photons travelling in a crossed-field 
configuration, where after we derive a new expression for the phase velocity of photons in a 
radiation gas. In the weak field limit, our new expression exactly reproduce the previously 
known results for the phase velocity, while the limit of an ultra-intense radiation gas gives 
rise to a constant subluminal phase velocity. It is argued that supcrluminal velocities are 
prohibited by pair creation. Our general dispersion relation shows that the effective force 
between the photons will always be attractive, similar to the weak field case. Thus, the 
collapse of inhomogeneities can occur within a perturbed gas of intense photons. Moreover, 
we investigate the implication for the collective dynamics of photons, using a kinetic theory. 
A set of equations governing the dynamics of an interacting gas and an electromagnetic pulse 
is derived. It is shown that they admit pulse self-compression and nonlinear saturation. The 
significance of our results are discussed, together with some possible applications. 

The effective field theory of soft photon-photon scattering in constant background elec- 
tromagnetic fields can be formulated using the Lagrangian density Jz? = Jzfo + -2c, where 
-£?o = —-^ is the classical free field Lagrangian, and the one-loop correction is [25] 

£> c = --t- / — e - m ' s {esfab coth(eas) cot(e6s) - \{es) 2 (a 2 -b 2 )-l , (1) 

where m e is the electron mass, e = |e| is the electron charge, a — [(,^ 2 + Sf 2 ) 1 / 2 + ^"J 1 / 2 , 
b = \{& 2 + Sf 2 ) 1 / 2 - J?] 1 / 2 , & ee F ab F ab /A = {B 2 - E 2 )/2, Sf = F ab F ab /A = -E B, 

pah = e abcd Fcd l 2 _ ThuS; g = ( fl2 _ b 2y 2 &nd |^| = flfo _ 

The dispersion relation for a test photon in a background electromagnetic field is given 
by [3] (see also Ref. [26]) 

[1 + \\{E 2 + B 2 )} v 2 - 2\k ■ (E x B)v + A \\{E 2 + B 2 ) - (k ■ E) 2 - (k ■ B) 2 ] - 1 = 0,(2) 

where v = u)/\k\ = 1/n is the normalized (by the speed of light) photon phase velocity, n is 
the refractive index, lu is the photon frequency, k is the photon wavevector, and k = k/\k\. 
The effective action charge A is [3] (see also [26]) 

\ = (d%+d 2 )J? 

-2dj?3? + &{d% + d 2 )^- 2(&d% + &d%^)3f ' ^ 1 

To start with, we look at crossed fields for which \E\ = \B\ and E ■ B = 0, i.e. a = b = 0. 
We define s = (E x B)/\E\ 2 , and assume that the test photon is a part of an ensemble of 
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random photons such that (k ■ s) = and ((k ■ s) ) = 1/3, where the angular bracket denotes 
an ensemble average. Then, Eq. J2J yields the average phase velocity v = [(1 — |A c f |-E| 2 )/ (1 + 
A c f l^l 2 )] 1 / 2 , where the index cf denotes the crossed-field value. Taking the limits a — ► b and 
b — > in Eq. and using the Lagrangian (JIJ, we obtain 

A -i = ^^ 4ritj (4) 

cf 22 Q cut, V I 

corresponding to held strengths ~ 8 x 10 17 V/cm or 3 x 10 15 G. Here, S CT n, = E% Tit> where 
E cr it = TO 2 /e = m 2 c 3 /eh ~ 10 16 V/cm denotes the critical Schwinger held and a — e 2 /An is 
the fine-structure constant. Note that this is the same charge as the geometrical average of the 
coefficient obtained from the polarisation tensor in the weak field limit [3], i.e. an average over 
polarisation states. Moreover, the Schwinger critical field will, in our case, not correspond to 
a threshold for the pair-creation. 

We will now analyse the radiation gas. An electromagnetic wave moving in an isotropic 
and homogeneous medium with the refractive index 1/v will satisfy \B\ — \E\/v, Sf = — E ■ B 
= 0, and — \{v~ 2 — 1)|-E| 2 > 0; the latter inequality holding for a subluminal phase 
velocity. This implies a ^ and 6 = 0. For supcrluminal velocities, the invariants satisfy 
a = and 6^0, for which the spontaneous pair production is predicted [25]. We will assume 
v < 1 throughout. A gas of photons may be viewed as an ensemble {Ei, Bi} interacting, due 
to elastic photon-photon scattering, through the refractive index of the gas. The effective 
action charge, as given by the expression will now depend on both the field strength and 
the phase velocity of the gas. But through the phase velocity itself depends on A, and 
is therefore implicitly determined via J2J). We assume that the momentum flux in the gas 
rest frame vanishes, i.e. (E x B) = 0, and the gas is characterized by its energy density 
$ = (|.E| 2 ). We, therefore, have a nonzero invariant according to a = [{v~ 2 — l)^] 1 / 2 . 

Equation J2J may now be solved to give the expression 



2 1 - f Agas^ + y 1 - 2A gas( f + i(A gas <?) 2 

2 + Agas^ 

for the phase velocity v. An expression of the same functional form can be derived for a random 
set of photons in the field of a plane wave pulse with S replaced by |-Ej,| 2 , where E p being the 
slowly varying pulse envelope. The field space Laplacian operating on the Lagrangian density 
yields (using z = aes) 

^9 % ^„ 1 e 2 f lco dz P _ /„ / 1 — z cothz 1 \ a 1 „, ,„ 
im {d%+dl)J? = — — / — e -E^z/a + thz = -F(a/E clit ), 



j-j5v-*--»/- 87r 2 2a 2 7 z \ sinh 2 z 3 J ~ 4tt a 2 ' 

(6) 

and gives the behavior of the effective action charge Q, while lim^o [— 29^-^f c — a 2 <9|?Jzf c ] = 
(a/8n)G(a/ E CI it) gives a small correction, and can safely be neglected. In its full form, Eq. 
(0 can be written as 



a F{a/E cHt ) 

W 2 + (a/8ir)[F(a/E CI i t ) + G(a/E c , it )} ' 



x _ « 1 yu/^ cl -itj . . 



The same functional relationship is found for a random set of photons in the field of a plane 
wave pulse, with § — > \E p \ 2 . The full normalized phase velocity v, as given by Eqs. ©-0, is 
plotted in Fig. as a function of the normalised energy $ /<? C rit- We see that the normalized 
phase velocity goes to unity for small energy densities, and decreases with increasing energy 
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Fig. 1 - The normalized phase velocity v, as given by Eqs. JSJ-Q, plotted as a function of the 
normalised energy density £ /£ cr it- 



densities. In the weak field limit {$ /# C rit <C 1), we take v ~ 1 in the evaluation of A gas , yielding 
Agas — A c f [see Eq. Q]. A series expansion of Eq. JSJ) yields v ~ 1 — §A c f#, consistent with 
the results in Refs. [3,14,26]. On the other hand, in the ultra-strong held limit (S j 1 §^\t 3> 1), 
F goes to infinity in Eq. J7J), giving A gas ~ 2a -2 = 2w 2 [(l — v 2 )S]~ x . Using Eq. (J5J, we obtain 
the asymptotic phase velocity 

Voc ~ 1/Vb w 0.45. (8) 

We stress that, formally, the one-loop approximation can be viewed as questionable in the 
limit § — > oo. However, it is straightforward to show, using higher order perturbation theory 
[27], that the renormalization group improved Lagrangian will include at most a logarithmic 
diverging behaviour, and thus the one-loop approximation gives a very accurate description 
even in the limit highly super-critical intensities (as in the case of pure magnetic fields [28]) 
Under most circumstances, the contribution proportional to a in the denominator of Eq. 
@ is small and may be neglected. When <S /<? C rit > 1, we obtain F ~ a/3E CI i t from Eq. JBJ, 
and we have 



A= 



8ira- 



:F(a/E„ it ) 



a 



1 



24-KE clit vT 



(9) 



Our knowledge of the nonlinear refractive properties of the radiation gas gives a means 
for investigating the effects of higher-order nonlinear corrections to the standard first order 
Heisenberg-Euler Lagrangian, and probe the significance of higher-order effects for photonic 
collapse [14, 15, 17]. The dynamics of coherent photons, travelling through an intense radiation 
gas, may be analysed following Ref. [14]. We obtain a nonlinear Schrodinger equation for the 
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slowly varying pulse envelope E p according to [29] 

| + v ■ V) E p + ^-V\E p + u Q ^fl Ep = o, (10) 
Ot J 2ko no 

where the subscript denotes the background state, V\ = V 2 — (feo • V) 2 , = $ — So is a 
perturbation, v — v(So), n = n(S ), n n \(8S) — J2m=i no8S m /m\, andr^" 1 - 1 = d m n /dS m . 

For a dispersion relation u> = \k\v(r,t), the motion of a single photon may be described 
by the Hamiltonian ray equations [30] 

r = dkUJ = v k, and k = — Vw = wV Inn. (11) 

Since n = n{S) and dnjdS > always holds, a denser region of the radiation gas will 
exercise an attractive force on the photon [31], thus creating lensing effects. The single particle 
dynamics thus supports that photonic self-compression is an inherent property of the one- 
loop radiation gas, but we note that as the density of a region increases, the phase velocity 
approaches a constant value JSJ), i.e. Vlnn-^0. 

Using Eq. (|11|) . the collective interaction of photons can then be formulated as [30] 

^ + vk.VN k + ^(Vn).^ = 0, (12) 
ot n ok 

where the distribution function Nk{r,t) has been normalised such that J Nk(r,t)dk is the 
number density. Interestingly enough, in the limit of large intensity S 3> fS'crit, since the phase 
velocity approaches a constant value, the one-loop radiation gas essentially evolves by free 
streaming, making higher-order loop corrections important in this extreme regime ( 1 ). 

Following Ref. [14], the response of the radiation gas to a plane wave pulse may be formu- 
lated in terms of an acoustic wave equation 

32 ^V 2 ) 5S = -^( ^ +v 2 vA n nl (\E p n (13) 



dt 2 3 J no \dt 2 

derived from Eq. 1|12[) by taking the moments, using the equation of state P = (a/3 (P being 
the pressure), and linearizing around the background state. If the time response of bS is slow, 
Eq. 1^31) gives 

+ (14) 
no \ 2n' J 

Using i|14|) and the expression for n n i((5cf), we can write Eq. (|10|) as 

2 



| + „ . v) E, + ^ L E r + m (2M) (l + Utf* = 0. ,15) 

When nQ\Ep\ 2 /2n' Q <C 1, we have a self- focusing nonlinearity in Eq. Ijl5|l . but as |-E P | grows 
the character of the nonlinear coefficient changes. The coefficient is positive when \E p \ 2 < 
Eg = |2no/no'|, but since n' ' < for all S'q, the sign will change as the pulse amplitude grows 
above the saturation field strength E$, making the nonlinearity defocusing and arresting the 
collapse. The numerical value of this turning point is dependent on the background parameter 
So- For low intensity radiation gases, n' ' ~ 0, and Eq. I|15|l always displays self-focusing, i.e. 

( 1 )The two-loop corrections are normally very small compared to the one-loop Lagrangian [27], making 
higher-order loop corrections important in the regime where v is close to the constant value J^J 
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the field strengths can reach values above the Schwinger field. When Sq roughly reaches the 
critical value <?crit, the weak field approximation breaks down, and Eqs. (J5J and 10 can be used 
to derive an expression for Eg- As an example displaying the general character of the intense 
background case, consider S = S a - lt X 10 2 . We find that Es ~ 2 x 10 17 V/cm > f^rlt, i.e. 
the pulse saturates above the Schwinger critical field. Thus, both the weak and moderately 
strong intensity cases, as described by Eq. I)15[i . display self-compression above the Schwinger 
critical field. The analysis of the saturation field strength for arbitrary background intensities 
requires that we abandon the approximation @ , and take the full denominator of Eq. J7J into 
account. 

An astrophysical setting where photon-photon scattering could result in photonic self- 
compression is the magnetar [19]. It has been shown that the quakes from crust tension, due 
to the extreme fields ~ 10 15 G, would release large amounts of low-frequency photons, creating 
a high intensity photon gas [32], suitably described within the framework developed here. The 
full consideration of such environments is beyond the scope of the present paper, and is left 
for future research. 

At the laboratory level, laser intensities are approaching 10 22 W/cm 2 and beyond [20]. 
While there are upper limits to powers achievable by ordinary laser techniques [20], laser- 
plasma systems hold the promise of surpassing these limits, coming closer to the Schwinger 
intensity 10 29 W/cm 2 [24]. Due to the relativistic ponderomotive force, the plasma particles 
will be expelled from the high intensity regions [33], and plasma channels will be created. As 
shown in Ref. [34], these channels provide the right conditions for photon-photon scattering 
to take place, and it has been suggested to use this a means for detecting elastic scattering 
between photons [35]. In the next generation laser-plasma systems, the radiation intensity 
within the plasma channels will reach levels where radiation vacuum collapse can take place. In 
such systems, the full nonlinear structure of the vacuum, as given by its refractive properties 
[Eq. ©], will be important. Moreover, it has been shown that once the pulse satisfies a 
collapse criteria, a self-compression of radiation gases will naturally occur due to QED effects 
[14, 15, 17, 18], making it possible to reach, within the self-consistent theory of elastic photon- 
photon scattering, the necessary energy scales for the novel strong-field effects to take place. 
As we have seen, the weak field approximation gives a good description up to Sq ~ <^ C rit, thus 
allowing self-compression of a plane wave pulse above the critical field strength -E7 c rit- O n the 
other hand, when £$ > <? cr it, the pulse will be saturated by higher order nonlinearities, but 
still above the critical field strength E a - lt . Thus, there is hope that one could probe the salient 
features of the nonlinear quantum vacuum in the near future. ( 2 ) 

To summarize, we have presented a new dispersion relation for a test photon in a radia- 
tion gas background, accurate to the one-loop level and valid for an arbitrary radiation gas 
intensity. We numerically analyzed the properties of the refractive index, and demonstrated 
photonic lensing due to energy density inhomogeneities. Furthermore, using the Hamilto- 
nian ray equations we developed a kinetic theory for photon-photon interactions, and found 
that the full nonlinear one-loop theory supports the notion of photonic self-compression up 
to a characteristic energy density where higher-order vacuum nonlinearities yield saturation. 
Moreover, the value of the saturated field strength was shown to depend on the radiation gas 
background energy density, and that it could surpass the Schwinger critical field strength. 

* * * 
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( 2 )It should be mentioned that derivative correction due to rapidly varying fields could be of importance in 
certain applications [36], e.g. involving 7-rays, but has little effect in the cases discussed here. 
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